Hysteresis loop signatures of phase transitions in a mean-field model of disordered 
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In accordance with recent experiments the mean-field type theories predict the presence of nu- 
merous metastable minima (states) in the rugged free-energy landscape of frustrated disordered 
magnets. This multiplicity of long-lived states with lifetimes greater than 10^ s makes the task to 
experimentally determine which of them has the lowest free energy (and thus what thermodynamic 
phase the sample is in) seem rather hopeless the more so as we do not know a protocol (such as 
field-cooling or zero-field-coohng) leading to the equihbrium state(s). Nevertheless here we show 
in the framework of Landau-type phenomenological model that signatures of the mean-field equi- 
hbrium phase transitions in such highly nonequihbrium systems may be found in the evolution of 
the hysteresis loop form. Thus the sequence of transitions from spin-glass to mixed phase and to 
ferromagnetic one results in the changes from inclined hysteresis loop to that with the developing 
vertical sides and to one with the perfectly vertical sides. Such relation between loop form and 
the location of global minimum may hold beyond the mean-field approximation and can be useful 
in the real experiments and Monte-Carlo simulations of the problems involving rugged potential 
landscape. Also the very existence of the quasi-static loops in spin glass and mixed phases implies 
that the known disorder-smoothing of the first-order transition can be always accompanied by the 
emergence of multiple metastable states. 

PACS numbers: 75.10.Nr, 75.50.Lk 



Many types of disorder can induce the frustration in 
the magnetic interactions of a solid. This results in ap- 
pearance of the glassy states in disordered magnets which 
may exist along with the more usual ones - ferromagnetic, 
ferrimagnetic etc. There are also mixed phases where 
perfect magnetic order is destroyed only partially. This 
variety of thermodynamic phases (states) is present in 
many microscopic models with frozen frustrated disorder 
[ij-Q. Yet an equilibrium thermodynamic phase may be 
something like a phantom in such highly nonequihbrium 
systems. From Monte-Carlo studies of microscopic spin- 
glass models Q it is known what a huge time is needed for 
such systems to relax into the equilibrium state. From 
the mean-field point of view this is the consequence of 
the rugged free energy landscape with numerous local 
metastable minima in each of which the system can be 
trapped on its way to the global (equilibrium) minimum 
for a time which grow exponentially with the volume of a 
sample J.] . Recent experiments d, Q show visually that 
in zero and finite field random magnets do stay confined 
in the variety of the phase space regions with different 
magnetizations for times up to lO^s without a noticeable 
relaxation. 

Moreover, Monte-Carlo simulations of short-range 3d 
Ising spin-glass under quasi-static field variations [l^l 
demonstrate the existence of numerous stationary states 
inside the hysteresis loop. These findings are shown to 
be in a reasonable qualitative agreement with the simi- 
lar experiments in the Ising spin-glass Fen^^IriQ^^TiO^ 
[loj . Along with the experiments of Refs.jsiQ this is the 
unambiguous evidence that short-range Ising spin-glasses 



do can posses the multi-valley structure of phase space 
even in external field. 

Often the validity of such mean-field picture for the 
short-range spin-glasses is questioned on the grounds of 
experiments (Ref. [ll| and references therein) and simu- 
lations 12|, [l3| which do not find any Almeida- Thouless 
anomalies in a field. It seems that the imaginary para- 
dox between evidences of the multi- valley structured - lOl 
and the absence of Almeida-Thouless anomalies |lll4l3| 
stems from the unjustified conviction that in short-range 
systems the onset of spin-glass phase should be necessar- 
ily accompanied by the divergence of the spin-glass sus- 
ceptibility same as in the Sherrington-Kirkpatrick model. 

So the rather natural resolution of this paradox lies 
in the recognition that in the short-range spin-glasses 
the multi-valley landscape can emerge in a finite field 
without the Almeida-Thouless anomalies. The qualita- 
tive explanation of such behavior can be found in the 
framework of multi-mode condensation mechanism 
Hij l describing the possible origin of the spin-glass phase 
in short-range random Ising magnets. Indeed, in this 
mechanism the whole macroscopic set of sparse fractal 
modes become unstable at the same transition temper- 
ature in zero field while in finite field their instability 
points spread throughout a finite temperature interval, 
cf. Fig. 4 in Ref. laj. So the experiments and simula- 
tions in the field would not detect the sharp transition, 
instead the gradual emergence of more and more corru- 
gated landscape would take place in the larger tempera- 
ture intervals for higher fields. 

Thus we may conclude that the mean-field viewpoint 
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is applicable to the real random Ising magnets at least 
on the present days' laboratory times. The last reserva- 
tion is important as usually the mean-field theories over- 
estimate the barriers between local minima and, hence, 
the Arrhenius lifetimes of the metastable states. Thus 
the Monte-Carlo simulations of short-range 3d Ising spin- 
glass it], |18| indicate that there can be the directions in 
phase space where the barrier's heights stay finite. Hence 
the mean-field notion of the divergent metastable states' 
lifetimes can be violated at longer time scales, that is 
the system can relax to the global minimum after, say, 
months or years. 

Yet, being in error with respect to the barriers and 
lifetimes, some refined mean-field models may still give 
correct relative depths of metastable minima in the free- 
energy landscape. In such cases we would have a correct 
mean-field description of a system's true phase diagram 
which would agree with the exact results of the equilib- 
rium statistical mechanics valid at infinite times. Thus 
the mean-field approach may be useful in the locating 
the points of the true equilibrium transitions. 

So with the above reservations we may interpret 
the experimental data in the mean-field regime spread- 
ing up to the times of order lO^s with the notion of 
metastable states with extremely long and divergent life- 
times. Hence, we are faced with the problem to find 
among the variety of the observed stationary states the 
equilibrium one(s) whose properties are described by the 
equilibrium microscopic models. In the absence of some 
notion of a possible protocol (such as field-cooling or 
zero-field-cooling) leading to the equilibrium state(s) it 
is rather hopeless task if one is not prepared to wait for 
months and years in (probably vain) hope that the mean- 
field predictions fail at longer times and our random sys- 
tem still find the way to the equilibrium state. 

Nevertheless the definite signatures of equilibrium 
thermodynamic transitions in some glassy magnets may 
be found in quite simple experiments. They can be en- 
coded in the temperature evolution of their hysteresis 
loop form 15|, |l9|. This is quite natural as the pres- 
ence of hysteresis loop is the main and unambiguous 
manifestation of the irreversibility appearing in glassy 
phases. Due to the mean-field picture the origin of this 
irreversibility lies in the presence of multiple metastable 
states so one can easily conclude that their magnetiza- 
tion curves should lie in the loop interior. This apparent 
conclusion needs no further theoretical justification once 
the existence of metastable states is adopted. Hence this 
should be the necessary feature of any adequate mean- 
field model of glassy magnets. 

In spite of its simplicity the notion of the loop filled 
with the metastable magnetization curves allows to ex- 
plain qualitatively the origin of all quasi-static irre- 
versible phenomena in glassy magnets. Thus assuming 
that in the spin-glass phase the observed inclined hys- 
teresis loops are filled with the metastable magnetization 



curves one can easily explain how the different values of 
physic al q uantities emerge in various experimental pro- 
tocols [l6j. For example, it becomes quite obvious that 
the field-cooled process brings the system to the upper 
outline of the loop while the zero-field-cooled one makes 
it join its lower outline (lif . 

Moreover the inspection of hysteresis loop form can 
give some notion on the relative depths of the metastable 
minima thus making a link with the predictions of 
the equilibrium statistical mechanics. On this way the 
quite plausible assumption is that the deeper zero-field 
metastable minimum the larger coercive field of its mag- 
netization curve (at which the minimum vanishes de- 
stroyed by the field opposite to the zero-field magneti- 
zation). As the ends of the metastable magnetization 
curves constitute the outline of the major hysteresis loop 
one can deduce from its form which of the metastable 
state has the lowest zero-field free energy thus being the 
equilibrium one. Such heuristic considerations for the de- 
termination of phase sequence in glassy systems were first 
advanced in Ref.[3|. Here we show that in the frame- 
work of the Landau-type phenomenological model the 
equilibrium phase transitions do show up in the hystere- 
sis loop form and this heuristics is exact in spin-glass and 
ferromagnetic phases and approximate in the mixed one. 

We consider the Landau potential for the magnetiza- 
tions of the sparse fractal modes rm, i = l,...,iVo 
of the form 



F(m) 



-f — [ml* — h [ml 



No 

Here [m*^] — Nq^ ^ m^. It is the simplest model suit- 

1=1 

able for our aim as at d > it has three magnetic phases 
(spin-glass, mixed and ferromagnetic) possible for the 
Ising-type ferromagnet with frustrated nonmagnetic dis- 
order. In the multiple local minima of F (m) rrii acquire 
just two values m+ > and m_ < so all their va- 
riety is parameterized by the quasi-continuous variable 
< n < 1 defining the fraction of the positive rrii at a 
given metastable minimum. The existence of these min- 
ima is limited by condition < and the stability con- 
ditions 



m+ (n) /2 < — m_ (n) < 2m_|_ (n) , x (n) 
Here 



_ dm (n) 
" dh 



m (n) = nm+ (n) + (1 — n) m_ (n) 



> 



(1) 



is the net magnetization of a sample in the metastable 
state and x ("■) is its magnetic susceptibility. Among 
these local minima those with the lowest free-energy 



3 



(equilibrium states) are defined by the conditions 



-0.5 (SG-M) and 



+ max ^— Tg, m (n)^^ + Sdm {nf' > 0, m (n) h > Q. 



The magnetization of equilibrium states nieq obeys the 
equation 



[t/ + max (-Tg, m^g)] TOe^ + dm\^ = /i 
and the parameter n — rieq of these states is 



correspond to Tg = — ro/(c — 1) 
Tg = -To/(c - d - 1) = -1 (M-FM). 
+ i^) + (") = 0, (2) jj^ SQ phase we have inclined loop ubiquitous in ran- 
dom magnets (FigdJ while in M phase the region of sta- 
bility of metastable states acquires the S-shaped form 
which results in the development of vertical jumps in the 
major loop (FiglJ]). This takes place right at SG-M equi- 
librium transition point and the critical fields ±hs and 
magnetizations ±ms at which the jumps start are defined 
through 



(3) 



1 H- sign {rrieq) \ min ( — 1 



(4) 



Also for the equilibrium value of spin-glass Edwards- 
Anderson order parameter 

q {n) ~ [iTT''^] ^ [itt]^ = 71 (1 — n) [m-^- (n) — to_ ('^)]^ 

we get from Eqs.(IT]),@, 

Qeq min(-Tg - rnlg,0) 

In zero field at Tg < Tf , Tg < we have 

■"^ - 0, n^q = 1/2, = -Tg 



^eq 



SO here the spin-glass (SG) phase is realized. At /i -> +0 
and [1 + d)Tg < Tf < Tg < the mixed (M) phase with 
partially ordered moments exists with 




(5) 



[Tf - Tg (1 + d)] /d 



At last for h — > -1-0 and t/ < {l + d)Tg < we have 
ferromagnetic (FM) phase where 



^eq 



1 



''eq 



1, 



qeq 



0. 



In spite of continuous variation of m, 



eqi 



''eq 



and 



ieq 



throughout all phases the transitions between them are 
of the first-order. This is because the SG and FM phase 
stay always stable at SG-M and M-FM transition points 
correspondingly. 

Now we turn to the temperature evolution of the hys- 
teresis loop accompanying these transitions. As both 
Tg and Tf are linear functions of temperature in a spe- 
cific sample there is a linear relation between them 
Tf — cTg + To with c, To = const. Choosing c > d + 1 
and To > we get the system with the SG-M-FM 
phase sequence when Tg diminishes. Figs HE] show the 
subsequent evolution of the hysteresis loop filled with 
metastable curves through such phase sequence. Here 
c = 5, To = 2d = 2so the transition temperatures 



3d 



3/2 



dmi 



The magnetizations ±ms are strictly proportional to 
the zero field equilibrium magnetization rueq (cf . Eq. ([5|) . 
Such transitions in the loop form are found in experi- 




FIG. 1: (color online) Hysteresis loop formed by metastable 
magnetization curves in SG phase at Tg = —0.45. Dotted 
lines show the stability limits of metastable states, dashed 
line indicates the lowest (equilibrium) minima. 




FIG. 2: (color online) The same as in Fig.l for M phase at 
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FIG. 3: (color online) Hysteresis loop in F phase at Tg 
-1.25. 



and in various random spin models [23 - |23 | 



ments 

We may reasonably suppose that they are the manifes- 
tations of the equilibrium SG-M transition in them. 

We can also mention another specific feature of M 
phase, that is the minor loops with vertical sides which 
tops (bottoms) are composed of metastable curves with 
sufficiently large n{l-n), see FigH At the M-FM tran- 
sition point TTis reaches the stability point of the homo- 
geneous state (n = I) so the loop in F phase becomes the 
ordinary ferromagnetic one composed entirely of homo- 
geneous magnetic states (Figl3]). The metastable curves 
inside it are inaccessible in quasi-static regime yet they 
can be reached with special fast field changes. 

The dashed lines in Figs. [TJ[5]show the position of the 
deepest equilibrium minima at a given field thus compris- 
ing the results for the equilibrium magnetization curves 
which can be obtained in the (mean-field) models of the 
equilibrium statistical mechanics. Apparently they are 
unobservable in the real-time experiments. But if one 
manage to get the experimental data in ac field with pe- 
riod much greater than lO^s the shrinking of loops to the 
equilibrium curves in Figs. [Tl[2]may be observed. If this 
will be actually the case and at what superlong times this 
may occur is still the open question of modern theory. In 
terms of the equilibrium curves the main feature of the 
present model can be characterized as simultaneous qual- 
itative changes of these curves and the loop's outline, e. 
g. the jumps at the loop's sides appear simultaneously 
with the jump in the equilibrium curve at h = indi- 
cating the onset of the mixed phase. The real magnets 
and more realistic models may not have such strict cor- 
respondence yet on the basis of above heuristics we may 
expect in them qualitatively similar behavior. We may 
also note that in random anisotropic Heisenberg magnets 
the evolution of loops would depend on the field direc- 
tion and the variety of mixed and spin-glass phases can 
be detected from this dependence. 

To elucidate the physical content of the present phe- 



nomenology we note that it still holds for the ordinary 
ferromagnets with low concentration of impurities and 
defects. Here the metastable states are just the numerous 
multi-domain patterns stable in some intervals of applied 
field. Each pattern has its own magnetization curve along 
which the net magnetization varies due to the change of 
the domain's moments. At the ends of this curve the 
pattern looses its stability and transforms into another 
one via domain wall jumps. The evolution of the hys- 
teresis loop filled with these curves in such slightly dis- 
ordered ferromagnet will follow that depicted in FigslT][3] 
as with lowering temperature the gradual growth the do- 
main wall energy will make lower the potential of the 
states with fewer walls and larger net moment. With in- 
creasing disorder the domains in these patterns will be 
smeared but nevertheless their spin configurations still 
will be stable corresponding to the local free energy min- 
ima. Thus the metastable spin-glass states can be viewed 
as smeared domain patterns. In the absence of domain 
walls the transitions between them are realized through 
the avalanches of spin upturns on the sparse fractal set 
of sites belonging to the specific eigenmode of random 
exchange matrix 
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16| . Hence Barkhausen noise and 
acoustic emission can be also observed at the sides of the 
outer loops in glassy phases. 

Figs HIS] shed also a new light on the smearing of the 
first order transitions by a quenched disorder [25[. In- 
deed, in perfect ferromagnet the field sweeping from neg- 
ative to positive values results in the first order transition 
in which magnetization direction is reversed. So Figs, [l] 
[2] show what smearing effect the disorder can have on this 
transition. Here the important feature consists in the ap- 
pearance of multiple metastable states in the vicinity of 
smoothed transition. It is quite probable that this is 
inherent property of all smeared first order transitions. 
Thus multiple magnetization curves show up near the 
AFM-FM transition in Fe-Rh alloy 

Here we can verify do the deepest (equilibrium) zero- 



field state with rrieq > {rieq > 1/2) has the largest 
module of its negative coercive field, i.e. the field de- 
stroying the local minimum. From Figs. [Tl[3]we see that 
it is evidently true in SG and F phases. Designating the 
parameter n of the state with the largest module of neg- 
ative coercive field as ric we have ric = n^q = 1/2 in SG 
phase and Uc = n^q = 1 in F phase. In M phase the ric 
of the state with largest coercive field is 

J 1/2, Tg (1 + 0.037d) <Tf < Tg 

^ I + v^^) ' + ^) + 

This Uc does not coincide with Ueq in Eq. ([5]) yet their 
values stay rather close throughout M phase as Fig. |4] 
shows. So here the above heuristics works only approxi- 
mately. 

Nevertheless the peculiarities of the major hysteresis 
loop form in the present model strictly follow the equi- 
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librium transitions. The thermodynamic phases in such 
system can be identified through simple experiments in 
the quasi-static periodic field varying sufficiently slowly 
to follow the potential profile. It can take much less 
time than the relaxation into the deepest minimum as 
encircling major loop the system sweeps the boundary 
between the rough and smooth landscapes and the relax- 
ation times here would not be that enormous. Yet the 
strict relation between loop form and equilibrium tran- 
sitions can be lost in more realistic models with com- 
plicated potential landscapes. One may argue that pos- 
sible deviations from the above picture should not be 
large as present model gives quite adequate de scrip tion 
of the loop evolution in real glassy objects and 
the underling heuristic considerations can still be valid 
for the real systems. Thus it is worthwhile to study with 



sis loop form as have been already noted in Ref. [23 



0.75 




FIG. 4: (color online) Temperature dependence of Ueq (solid 
line) and ric (dashed line) in M phase 

Monte-Carlo simulations the relations between the loop 
form and the global minimum position in the microscopic 
short-range models. Moreover just these studies can save 
the results of the equilibrium statistical mechanics and 
Mont-Carlo absolute minima searches from being use- 
less for the real experiments. Figs[l][3] show unambigu- 
ously that full knowledge of the potential landscape and 
its dependence on the field and temperature are neces- 
sary for understanding the quasi-static processes in glassy 
systems. It would not be a formidable task for experi- 
ments and simulations as it just comes to the determina- 
tion of the temperature evolution of hysteresis loop and 
metastable curves inside it. In reward one gets the ability 
to describe every (unavoidably nonequilibrium) process 
on the experimental time scales [H, [3] . 

Also this strategy and heuristic relation between the 
coercive fields and minima depths can be used in the 
variety of optimization problems once a suitable couple 
of conjugate variables playing the role of m and h can 
be found. For example, it is quite evident from Figs. 
[U [5] that the ground state search with numerical ver- 
sion of demagnetization process f27| can converge to the 
non-magnetized ground state but should be modified for 
the magnetized one. So the sensible ground state search 
should include preliminary determination of the hystere- 
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